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1 Introduction 

1.1 Summary of results 

In this paper we study the connection between dimers and Harnack curves 
discovered in ^Hj- To any periodic edge- weighted planar bipartite graph F 
one associates its spectral curve P(z, w) = 0. The real polynomial P defining 
the spectral curve arises as the characteristic polynomial of the Kasteleyn 
operator in the dimer model. 

It was shown in ^3] that for real positive edge weights on F the curves 
thus obtained are real curves of a very special kind, namely they are Harnack 
curves. Harnack curves are, in some sense, the best possible real plane curves. 
They were studied both classically and recently, see [TBJ HI] and references 
therein. Here we prove that every Harnack curve arises as a spectral curve 
of some dimer model. This gives a parameterization of the set of Harnack 
curves which, in spirit, is similar to the classical parametrization of totally 
positive matrices. It may be compared with the result of Vinnikov [21J who 
gave a similar description of real plane curves with a maximally nested set 
of ovals (which form a class of curves in some sense opposite to Harnack 
curves) . 

We also prove that modulo the natural (R x ) 2 -action the set of degree d 
Harnack curves in P 2 (IR) is diffeomorphic to the closed octant ]R>q ^ d 1 . 
In fact, the areas of the amoeba holes and the distances between the amoeba 
tentacles give these global coordinates. 

The Kasteleyn operator of the dimer model is an example of a periodic 
finite-difference operator (weighted adjacency matrix of a periodic graph). 
The spectral theory of such operators is much developed, stimulated, in 
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particular, by their connections with integrable systems, see, for example, 
[21 El M f° r an introduction. The particular type of operators we consider 
was studied by A. Oblomkov ^H] and in a series of paper by I. Dynnikov and 
S. Novikov pHE]. The spectral data associated to a periodic finite-difference 
operator is, typically, an algebraic variety (a curve C, in our case) together 
with a line bundle on it, that is, a together with a point of the Jacobian 
J(C). With complex coefficients, it is usually easy to see (see e.g. [5]) that 
the "spectral transform" , from the operator to its spectral data, is dominat- 
ing, but not surjective. Reality issues are typically more subtle. From the 
probabilistic origin of our spectral problem, it is natural to consider real and 
positive coefficients (edge weights). This adds a new aspect to the problem. 

The Harnack curves of genus zero play a special role. We characterize 
them as the spectral curves of isoradial dimers studied in ^3], see Section 
and also as those Harnack curves that minimize the volume under their 
Ronkin function (with given boundary conditions), see Proposition El Trans- 
lated into the language of probability, this means that isoradial dimer weights 
maximize the partition function with given boundary conditions. 

1.2 Acknowledgments 

This research was started at the Institut Henri Poincare and completed while 
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DMS-0096246 and a fellowship from Packard foundation. 

2 Background 

2.1 Kasteleyn operator and its spectral curve 
2.1.1 

Let T be a periodic planar bipartite graph. We can assume that V is embed- 
ded in the plane M? in such a way that translations by the standard lattice Z 2 
preserve T, including the partition of its vertices into black and white ones. 
Throughout this paper, we assume that the quotient T/Z 2 is finite. Let the 
edges of T be weighted in a Z 2 -invariant way. 

A perfect matching (also known as a dimer configuration) of a graph is a 
collection of edges that cover every vertex exactly once. The weight of such a 
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matching is denned to be the product of the edge weights. Kasteleyn in [T2] 
computed the weighted number of perfect matchings of any finite bipartite 
planar graph Tq using the determinant of what is now called the Kasteleyn 
operator 

K ' (£jklack vertices ^ ^ white vertices -j\ 

which is the weighted adjacency matrix of T sign-twisted in a certain way, 
see [T3j for an introduction. 

The study of perfect matchings (also known as the dimer model) on in- 
finite periodic planar graphs T stimulates the study of the corresponding 
Kasteleyn operators, see for example for further information. In this pa- 
per, we will focus on the case when T is the hexagonal lattice, in which case 
K is simply the weighted adjacency matrix. This case is, in fact, universal as 
will be explained in Section 12.1.41 



2.1.2 

The operator K commutes with the translation action of Z 2 and, in particular, 
it preserves the Z 2 -eigenspaces. These eigenspaces are indexed by characters 
of Z 2 , that is, by a pair of Bloch-Floquet multipliers {z,w) G (C*) 2 . They 
are finite-dimensional with a distinguished basis {5 V } consisting of functions 
supported on a single Z 2 -orbit and taking value 1 on a vertex v inside a fixed 
fundamental domain. Let K(z,w) be the matrix of K in the basis {S v } and 
acting on the (z, w)-eigenspace of Z 2 and set, by definition, 

P{z,w) =detK{z,w) . (2.2) 

Different choices of the fundamental domain lead to polynomials that differ 
by a factor of the form z l wK In particular, the zero locus of P 

{P(z,w) =0} C (C*) 2 (2.3) 

is defined canonically and is called the spectral curve of the dimer problem. 

The spectral curve remains the same if we multiply the weights of all edges 
incident to a given vertex by the same number. This is a gauge transforma- 
tions of the dimer problem; it does not change probabilities of configurations. 

The map from the weights modulo gauge to the corresponding spectral 
curves is our main object of study in this paper. 
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2.1.3 



For example, for the hexagonal lattice with (ixd-fundamental domain, P(z, w) 
is the determinant of the d 2 x d 2 matrix the construction of which is illus- 
trated in Figure [T] We denote the weights of the 3 edges incident to a white 
vertex v by a v , b v , and c v . The dashed line in Figured] is the boundary of 
the fundamental domain. An edge crossing it is weighted by an extra factor 
of z or w. 






b v 








Figure 1: The operator K(z,w) 



It is clear that in this case P(z,w) is polynomial in z and w of degree 
at most d in each variable. In fact P has total degree d: this can be seen 
by splitting the fundamental domain in half into two equilateral triangular 
arrays of vertices; the upper right triangle has d more white vertices than 
black vertices so exactly d edges connect this part with the rest of the graph. 

2.1.4 

The case of the hexagonal lattice is universal in the following sense. By 
choosing d large enough and setting some edge weights to zero, one can 
produce a graph which is equivalent to an arbitrary planar periodic bipartite 
graph, in the following sense. 

We consider two graphs to be equivalent if one can be obtained from 
another by a sequence of moves of the following type: remove a 1-valent 
vertex and its neighbor, or remove a 2-valent vertex, gluing its neighbors 
into a single vertex and redistributing the edge weights accordingly. See ^H] ■ 
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It is easy to see that such transformation induce a weight-preserving bijection 
of sets of dimer configurations. 

Because of this universality we will focus in this paper on the case of 
hexagonal lattice with d x ci-fundamental domain. In this case, the spectral 
curve is a degree d curve in the projective plane P 2 . 

2.2 Harnack curves 

2.2.1 Topological configuration of ovals 

A real algebraic curve C(R) C P 2 (R) is called an M-curve if it has the max- 
imal possible number of connected components, namely 1 + (d — l)(d — 2)/2, 
where d is the degree of C. For brevity, we will call all these components 
ovals, even though this is at odds with the classical distinction between sep- 
arating and nonseparating components of C(R) (which will play no role in 
this paper). We will also treat isolated real points as degenerate ovals. Ovals 
not intersecting the coordinate axes will be called compact ovals. 

Among the M-curves, the Harnack curves are distinguished by the posi- 
tion of their ovals with respect to the each other and the three coordinate 
lines in P 2 (they are also sometimes called simple Harnack curves in the liter- 
ature). The classical definition of a Harnack curve is somewhat complicated, 
see [123 . We will use instead modern characterizations of Harnack curves 
obtained in [TB*1 117] • By the main result of [16] . for given degree d, a curve 
is a Harnack curve if and only if the topological configuration of ovals is 
as illustrated in Figure El (It is different for odd and even degree.) The 
numbers of ovals in each quadrant are consecutive triangular numbers. Of 
course, Figure 121 is only meant to illustrate the topology and not be a plot of 
an actual Harnack curve. 

The group (R x ) 2 acts freely on the set of Harnack curves by rescaling 
the variables. We will call the quotient of this action the moduli space of 
Harnack curves. We will see that the moduli space of degree d Harnack 
curves is diffeomorphic to the closed octant M^q^ 1 ^ 2 , in particular, it is 
connected and simply connected. It follows that the set of all Harnack curves 
has 4 connected components, corresponding to a choice of one of 4 quadrants 
in Figure 121 These choices are related by reversing the signs of z or w. 

There is a corresponding flexibility on the dimer problem side: rescaling z 
and w corresponds to a natural operation on dimer weights which was called 
changing magnetic field in [TSj. From many points of view, it is natural 
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Figure 2: Oval configuration of Harnack curves of degrees 6 and 5 



to view this operation as a generalized gauge transformation. For example, 
with such extended definition of gauge equivalence, any dimer model whose 
spectral curve is rational is equivalent to an isoradial dimer, see Section 

2.2.2 Amoeba map 

Another useful characterization of Harnack curves C is that the map 

(C*) 2 D C(C) 3 (z,w) i-> (log |4 log \w\) e A{C) c R 2 

from the curve C to its amoeba A{C) is generically 2-to-l over the interior 
of A(C) (that is, 2-to-l except at real nodes). 

The (geometric) genus g of a Harnack curve is the number of compact 
ovals that are not reduced to points. In particular, the amoeba of a genus g 
Harnack curve has exactly g compact holes, which is illustrated in Figure 03 
for the case g = 3. 

The 2-to-l property constrains the possible singularities that C can have 
by constraining the possible links of the singular points. It was shown in j^j 
that the only possible degenerations of Harnack curves with fixed Newton 
polygon occur when some of the ovals shrink to zero size, producing real 
isolated double points (real nodes). In particular, it is impossible for two 
ovals of C to meet. 
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2.2.3 Ronkin function and Monge- Ampere equation 

A third characterization of Harnack curves C, obtained in is by maxi- 
mally of the area of their amoeba A(C). Concretely, C is Harnack if and 
only if its amoeba has the maximal possible area for given Newton polygon 
A, namely 7r 2 times the area of A. 

This maximality of the area is a consequence of the following Monge- 
Ampere equation (see [T7| ) 

det (R xx M=i_, (2.4) 



>*yx ± ^yy 



satisfied for all (x, y) in the interior of the amoeba A(C). Here R(x, y) is the 
Ronkin function of C defined by 

R( X} y) = 7 J— [[\og\P(z,w)\ -—■ (2-5) 



(27Ti) 2 / J\z\=e x ' z w 

y ' • J - J \ w \= e y 
2.2.4 Holomorphic differentials 

Holomorphic differentials on C will play an important role in our analysis. 
In this section we recall some basic facts about them. 

Suppose that C has geometric genus g. Let a%, . . . , a g be the compact real 
ovals of C that are not reduced to points, and let {g^}, i — 1, . . . , ("^j — 9 
denote its isolated real nodes. The ovals aij form a set of a-cycles in Hi(C). 
The corresponding 6-cycles can be taken in the form 

/3 i = {\z\ = e Xi ,\w\<e Vi }nC (2.6) 

where (xi,yi) is a point inside the ith hole of the amoeba A(C). The cycles 
Pi are anti-invariant with the respect to the complex conjugation. 
Holomorphic differentials on C have the form 

Q(z,w) 

where Q is a polynomial of degree d— 3 vanishing at the points see pQ. 
It is a fundamental fact about plane curves (see pQ, Appendix A) that the 
nodes {qi} impose independent conditions on polynomials of degree d — 3. 
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In particular, the space of differentials (|2.7|) is g-dimensional. We choose a 
basis {uji} such that 

/ V (2-8) 

The polynomials Qi are real and their 6-periods are purely imaginary (since 
Pi is anti- invariant under complex conjugation). 

The polynomial Qi has at least two zeros on every compact oval except 
the ith, by the integral condition ()2.8|) : the form dz/P w is real and of constant 
sign on every oval, so Qi must change sign at least twice on each oval except 
the zth. This implies that Qi must have no zeros on the i-th oval and precisely 
2 zeros on any other oval, otherwise, the degree d — 3 curve Qi = would 
intersect C at 

(d-l)(d-2) > d(d-3) 
points. It follows that the function 

ut, (2.9) 

where Xq and x are points on the zth oval, is increasing and provides a natural 
parametrization of the ith oval. 

Let us call a degree ( <l! 2 1 ) divisor D a standard divisor if it has precisely 
one point on each compact oval of C (including degenerate ovals). By mono- 
tonicity of (|2.9|) . the product of compact ovals of C, which is the variety 
parameterizing standard divisors, injects into the Jacobian J(C) of degree 
. It forms a component of its real locus. 



3 The spectral problem 

3.1 Boundary of the spectral curve 

By the boundary of a spectral curve C we mean the 3d points (counting 
multiplicity) of its intersection with the coordinate lines of P 2 . We show 
in this section that these points have a simple interpretation in terms of 
the weights of the dimer problem. In particular, they are real for real edge 
weights. 

It is useful to introduce the curve C defined by the equation 

P(z,w) = P(z d ,w d ) . (3.1) 
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This is a <i 2 -fold covering of the spectral curve C ramified over its boundary. 
The advantage of going to this branched covering is that operator K(z d ,w d ) 
is gauge equivalent to the operator K(z,w) defined by the rule illustrated 
in Figure H3 In the operator K(z,w), the weights of all vertical edges are 
multiplied by w and weights of all northwest-southeast edges are multiplied 
by z. 




Figure 3: The operator K(z,w) 



In particular, if w = then this operator becomes block-diagonal with 
blocks corresponding to chains of vertices of the form shown in Figure |U 



/Y Y T / 



Figure 4: The operator K(z, 0) 
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The determinant of the corresponding block equals 



det 



\za d 



za 2 

c 2 za 3 
c 3 



J 



(3.2) 



Therefore, the spectral curve P(z, w) 
points of the form 



intersects the w = line at the 



:-ir 



n 



(3.3) 



v in a zig-zag cycle 



where by a zig-zag cycle we mean a "horizontal" chain of vertices like in 
Figure |U 

The alternating products of the from ()3.3|) taken for all paths of all 3 
orientations multiply to 1. This is, in fact, a general constraint for the 
boundary points of any degree d curve. It can be viewed, for example, as 
coming from the divisor class group of the union of the coordinate lines. 

Note that once the fundamental domain is fixed, we get a canonical order- 
ing of the boundary points of the spectral curve. It comes from the ordering 
of the corresponding zig-zag cycles. 



3.2 The divisor of a vertex 

By definition, the spectral curve ()2.3|) is where the kernel and the cokernel 
of the matrix K(z, w) become nontrivial. For smooth spectral curves C, the 
kernel and the cokernel are 1-dimensional [5 j and form a line bundle over C. 
In particular, for every white vertex v G T the image of the function 5 V in 

Coker K(z, w) = C white vertices / Img K(z, w) 

defines a section of the cokernel bundle. We will denote by (v) the divisor of 
this section excluding the points on the boundary. 

In down to earth terms, equations of (v ) are the cofactors of the entries 
of K(z, w) in the row corresponding to the white vertex v. Since the columns 
of K(z,w) are linearly dependent on the curve C, (v) is defined locally by 
the vanishing of a single minor. 
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Theorem 1. For real positive edge weights, the spectral curve P(z,w) = is 
a Harnack curve of degree d, in particular, it has ( n 1 ) compact ovals. The 
divisor (v) of any vertex v is a standard divisor on this curve. 

The first part of this theorem is one of the main results of [T^]. The 
argument given below provides a new proof of it. 

Proof. It suffices to consider the case of generic edge weights, since the prop- 
erties of being a Harnack curve and being a standard divisor are closed. 

Let us consider the bundle £ formed by the cokernels of the matrix K(z, w) 
on the curve C defined by the equation (|3-H) . It fits into the following exact 
sequence 

-> 0(-lf d2 -> £ -> (3.4) 
of sheaves on P 2 . It follows that 

The curve C has genus (d 2 — l){d 2 — 2)/2 and, therefore, by Riemann-Roch 
we obtain 

d\d 2 - 1) 

deg L = . 

& 2 

Let v be a white vertex and consider the section s v of £> defined by the 
function 5 V . The divisor of s v is the pull-back of (v) except for the boundary 
contribution, which we will now determine. 

Let p be a boundary point of C corresponding to a zig-zag cycle as in 
Figure Eland suppose that the vertex v is k steps above this path. Below we 
prove the following: 

Lemma 1. The order of the vanishing of s v at the point p equals k. 

Since every path gives rise to d points on the boundary of C, it follows 
that the total contribution of the boundary points to the degree of s v equals 
3d 2 (d — l)/2. Subtracting it, we obtain 

1 / n , \ (d-l)(d-2) , N 

deg(v) = ^(deg£ - degs v ) = K - } j >- . (3.5) 

The remainder of the proof is based on the deformation to the constant 
weight case. If all edge weights are equal to one, the polynomial P takes a 
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particularly simple form (see |15| ) 



^[z, w) = Yl ( £i z + £j w + !) > ( 3 - 6 ) 

where e is a primitive <ith root of unity. In particular, the spectral curve C 
is a genus curve with ( d ^ > 1 ) isolated real nodes (isolated real solutions to 
e l z + w + 1 = exist when i ^ j and i,j ^ d; moreover replacing i,j 
with gives an identical solution). It also intersects each coordinate 

line once with multiplicity d, and, hence, satisfies the topological definition 
of a Harnack curve. 
The space 

Ker K*{z,w) = (Coker K{z, w))* 

is 2-dimensional if (z, w) is one of the nodes of C. In particular, there exists 
an nonzero element of this space which annihilates 5 V . By continuity, for all 
nearby curves the divisor (v) has a nearby point, possibly complex. By (|3.5|) 
the total number of such points equals the degree of (v), hence all of them 
are simple real zeros (no point can bifurcate into a complex conjugate pair 
of zeros). 

It follows that under a small real perturbation of weights, the spectral 
curve remains a Harnack curve and (v) remains a standard divisor. Because 
the ovals of a Harnack curve cannot meet (see Section l2.2.2|) and the boundary 
points do not become complex (see Section l3~Tj) . the same statement holds 
globally for real positive weights. 

□ 

Proof of Lemma{]\ Let (z ,0) be the coordinates of the point p. Note that 
if the corresponding intersection of the curve C with the w = axes was 
transverse (which happens generically) then w is a local parameter on C 
near p and z — z = 0(w d ). Therefore, near p we have 



K 



c 



( K± wB 1 \ 

K 2 wB 2 



\wB d K d ) 



+ 0(w d ) 



where the blocks Ki correspond to zig-zag paths as in ()3.2j) and Bi are diag- 
onal invertible matrices. The matrix K\ has a one-dimensional kernel while 
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all other K^s are invertible. Clearly, 

Ker K* = (v, -wK^ l B lV , w 2 K^ 1 B 2 K^ l B 1 v, . . . ) + 0(w d ) , 

where v G Ker K*. Since the entries of Ker K* are the cofactors of K that we 
need, the lemma follows. □ 



3.3 Spectral transform 

We have constructed the following correspondence 

/ Harnack curve \ 

ordering of boundary points 



/ edge weights / gauge \ 

fundamental domain 
V vertex / 



(3.7) 



y standard divisor J 



which we will call the spectral transform. Our first result about it is the 
following 

Theorem 2. The spectral transform is injective. 

Proof. The spectral data determines the sheaf £ introduced in the proof of 
Theorem n up to isomorphism. By Theorem 1.1 in [5], this determines the 
matrix K(z,w) up to multiplication on the left and on the right by elements 
of GL(d 2 , M), that is, up to a choice of linear basis of both spaces in (|2.1j) . 

Let p be a boundary point of the spectral curve. It determines a filtration 
on the space C whlte vertlces by the order of the vanishing of the corresponding 
section of the cokernel bundle at p. It is easy to see that the basis {6 V } is 
the unique, up to normalization, basis compatible with all these nitrations. 
Same argument applied to transpose matrix K*, which has the same spectral 
curve, reconstructs the delta-function basis of C black vertlces . As a result we 
reconstruct the matrix K(z,w) up to a multiplication on the left and on the 
right by a diagonal matrix, that is, up to a gauge transformation. □ 

A very useful (and well-known, cf. [SI HDl HH |20j) way of thinking about 
the spectral data in ()3.7|) is the following. The divisor is an ordered collec- 
tion of distinct points of P 2 ; similarly, the boundary points are an ordered 
collection of points on the coordinate axes. The spectral curve P(z,w) is 
determined by this data as the unique degree d curve passing through all 
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these points. Indeed, since the boundary points are fixed, any other curve 
passing through the same points is given by 

P(z, w) + zwQ(z, w) = , deg Q < d — 3 . 

But then, as in Sect ion |2.2.4[ Q = and P = will have (d — l)(d — 2) points 
in common, which is one too many. 

It is clear from this description that the spectral data is parametrized by 
an open subset of R d2+1 (which, in fact, is homeomorphic to R d2+1 , as we will 
see later). Also, this picture continues to work in a complex neighborhood 
and the injectivity remains valid with the same proof. It follows at once that: 

Corollary 2. The differential of the spectral transform is injective. 

3.4 Surjectivity of the spectral transform 

This reconstruction procedure given in Theorem |2 can be made explicit to a 
certain degree, but it appears to be difficult to establish the surjectivity of 
the spectral transform directly. Instead, we will take a different route based 
on the following 

Proposition 3. The spectral transform is proper. 

In fact, at several points in the paper we will resort to proving surjectivity 
by combining properness with injectivity of the differential. 

Proof. Suppose that the spectral curve C varies in some bounded set of 
Harnack curves. In particular, the the coefficients of 1, z d , and w d in the 
equation of C are bounded from above and below in absolute value. Without 
loss of generality, we can assume that these coefficients are equal to ±1. This 
means that the weights of the 3 frozen matchings (matchings using all edges 
of the same orientation) are equal to 1. Other coefficients of P(z,w), which 
enumerate other periodic matchings, are bounded from above in absolute 
value. This implies that the weight of any periodic matching is bounded 
from above. 

Periodic matchings are the vertices of the periodic flow polytope. This 
polytope is formed by nonnegative periodic flows from white vertices to black 
vertices such that the total flux is 1 at each black vertex. The logarithm of 
the weight of a matching extends as a bounded linear function to the interior 
of this polytope (recall that the weight of a matching is the product of its 
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edge weights). The barycenter of this polytope is the flow of intensity 1/3 
along each edge. It has weight 1 by our assumption. 

Coordinates on the space of weights modulo gauge are given by the com- 
binations of the form 

weight (e 2 i-i) 



n 



weight (e 2i 



(3.8) 



where ej, e^, ■ ■ ■ , t2k is a closed loop of edges on the torus, such as for example 
the loop in Figure |U Such loops can be thought of as flows with zero flux. 
If the weight (13. 8|) of any such loop under the variation is unbounded, we 
can add a small multiple of it to the constant 1/3 flow and obtain a point 
inside the periodic flow polytope with unbounded weight. This contradiction 
completes the proof. □ 

We will prove in Theorem 03 that the spectral data on the left in ()3.7j) 
forms a manifold diffeomorphic to IR d2+1 . From the description (J3.8|) of the 
coordinates on the space of weights modulo gauge, it it clear that this is 
also diffeomorphic to M d2+1 . Indeed, the cycle space is d 2 + 1-dimensional, 
generated by the d 2 faces in a fundamental domain (subject to one relation), 
and the horizontal and vertical cycles. Therefore, we obtain 

Theorem 3. The spectral transform is a bijection. 



4 The space of Harnack curves 
4.1 Harnack curves of genus zero 

Harnack curves of genus zero can be understood explicitly. By general theory, 
any real degree d genus zero curve C has a parametrization of the form 

z(t) = a f[ t -^-, ^ = 4^' 

i=l i=l 

where t £ KLP 1 = iU {oo}, a$ and bo are nonzero real numbers, and the 
other parameters are either real or appear in complex conjugate pairs. Such 
representation is unique up to the action of PGL 2 (M) by fractional linear 
transformation of the variable t. 

The action of PGL2(M.) can by reduced down to the action of the con- 
nected group PSLi2(M.) once we fix the orientation of the curve C. We will fix 
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some cyclic ordering of the coordinate axes and require that the C intersects 
them in that order (see Section |2.2.1|) . Let the cyclic order be: {z = 0}, 
{w = 0}, and the line at infinity. 

Proposition 4. The curve (j4.1|) is Harnack if and only if the parameters in 
(|4.1|) are real and cyclically ordered on WP 1 as follows: 

a\ < a 2 < • ■ ■ < ad < bi < • ■ ■ < bd < c\ < ■ ■ ■ < Cd < ai . (4.2) 

Proof. In one direction, this is a part of the topological definition of the 
Harnack curve, see which requires the curve C to first intersect the line 
{z = 0} d times, then the line {w = 0} also d times, and, finally, the line at 
infinity d times. 

In the other direction, we have to show that as long as the parameters 
(14. 2|) are cyclically ordered the curve C remains Harnack. Note that the strict 
inequalities in ()4.2|) ensure that C doesn't pass through the intersections of 
the coordinate axes, which in other words means that the Newton polygon 
of its equation remains the full triangle. Since the real nodes of the curve 
cannot disappear or merge by the 2-to-l property (see Section 12. 2. 2|) . and the 
region in the parameter space defined by the inequalities (j4.2j) is connected, 
the result follows. □ 

The PSX2(IR)-quotient can be taken by fixing, for example, 

ai = 0, b\ — 1, ci = oo . (4.3) 

This gives the following: 

Corollary 5. The moduli space of genus degree d Harnack curves is dif- 
feomorphic to M> d _3 . 

Let us associate to the curve C its 3c? points of intersection with the 
coordinate axes, counting multiplicity. For the curve (j4.1)l these are the 
points in projective coordinates (0, Ai, 1), (1, 0, Bi), and (Ci, 1, 0) where 

4 = &dII — > ^ = 1 II^' c > = TU c -^r- ( 4 - 4 ) 

- LJ - ai — Cj an bi — a^ On Q — bj 

J 3 3 

Note that all Ai have the same sign (similarly for B^s and CVs) and the 
following relation: 

]jA i B i C i = (-l) d . (4.5) 
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Modulo the action of (M x ) 2 , these constraints define a manifold with bound- 
ary diffeomorphic to M> d -3 . To see this, sort the AJs and look at the ratios 
between consecutive values. This d— 1-tuple of , A' ratios is in M^ 1 . Similarly 
for the and 'C ratios. For a given set of ratios the values of the smallest 
Ai,Bi, Ci are then determined up to (M x ) 2 by (J4.5j) . If on the other hand we 
choose an ordering of the boundary points as in (J3.7j) then the constraints 
define an M 3d ~ 3 : for example by the (M x ) 2 -action we can fix A\ = B\ = 1 
then all other Aj, Bi, Ci except C\ are free, G\ being fixed by (|4.5|) . 
We have the following: 

Theorem 4. There exists a unique genus zero Harnack curve for every choice 
of boundary points. 

Proof. The map (J4.4)) descends to a map 

K bi, Cl } l=1 ... d /PSL 2 (R) -> {A t , B t , C t } t=1 ... d /(R X ) 2 , (4.6) 

which becomes a map from M 3d ~ 3 to IR 3d ™ 3 if we introduce an ordering of 
points on both sides. It is easy to see that this map is proper. For example, 
in the gauge ([4.3)1 . we clearly have bj — Cj > 1 for any j and any % > 2. If the 
ratio 

Ci t— r CLj Ci 

Cl ~ Al ^ - a 

is bounded from below then a,j — Ci is bounded from below. The same argu- 
ment (using the fact that aj — Ci and hence — q is bounded below) shows that 
if Ci/Ci and A4/A1 are bounded below then bj — a>i is bounded from below. 
Finally if Bi/ Bi is bounded from below then 



B\ f *- bi — aj f A 1 — c,- bi 

J I j=2 



and we conclude that the b^s are bounded from above. 

It suffices to check, therefore, that the differential of the map (J4.6)) is an 
isomorphism at every point. We compute 

dlogA, W 1_ 1 



j 



dak \o-i — bj Cj — ai 



gl °S^ 1 _ dlogB^ e ^ 

96j aj — 6j 9ai 
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Note that the Jacobi matrix (J4.7)) is a symmetric matrix which is the sum of 
the following elementary 3x3 blocks: 



a>i - bj 



+ 



1 



CI; 



1 



b, 



+ 



Ck - a-i 
1 



\ 



di - bj bj - c k 
1 



1 



b, 



Ck 



\ 



bj - c k 



+ 



1 



(4.8) 



Ck - CLi 

over all triples of indices i,j,k = l,...,d 
with kernel spanned by the vectors (1, 1, 1) 
eigenvalue equals 

(oi - bj) 2 + (c fc - cii) 2 



bj - c k ) 



Ck O'i 

The matrix (|4.8|) has rank 1 
and (04, bj,Ck). The remaining 



+ {bj - c k ) 



{cii - bj){c k - ai)(bj - c k ) 

which is nonvanishing and of the same sign for all triples k) by the cyclic 
ordering condition (|4.2|) . 

It follows that a vector is in the kernel of the Jacobi matrix (J4.7|) if and 
only if it is annihilated by every 3x3 matrix ()4.8|) . It follows that this kernel 
is spanned by (1,1,..., 1) and 

(ai, . . . , <2d,bi, . . . ,bd, ci, . . . , Cd) 

and coincides with the tangent space to the orbit of the subgroup 23 C 
PSL 2 (M.) formed by upper-triangular matrices. Indeed, replacing t i— >■ at + (3 
in (14. lj) corresponds to the change a; i— > (a^ — (3) / 'a and similarly for 6j and 

Ci- 

Since the matrix (|4.7j) is symmetric, its image is the orthogonal com- 
plement of its kernel. Now the tangent space to the PSL 2 (R)-orbit is 3- 
dimensional, contains the kernel of the differential of (|4.4J1 . and is mapped 
by this differential to the tangent space to the (M x ) 2 -orbit. We need to know 
the intersection of the image with the tangent space to the (M x ) 2 -orbit. It 
is immediate to see that this intersection is always 1-dimensional, in fact 
spanned by the vector 

(A, . . . , A, 5, . . . , 5, -A - 5, . . . , -A - 8) 

where \ Y^ a i + = (A + 5)J2 c i- It follows that this one-dimensional 

space has to be the image of the tangent space to the PSX 2 (K)-orbit. This 
shows that the differential of (J4.6|) is an isomorphism, which concludes the 
proof. □ 
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4.2 Intercept coordinates 



Let C be a Harnack curve with the equation P(z, w) = 0. The Ronkin func- 
tion (|2.5|) of the polynomial P has a facet (that is, a domain on which is it 
affine linear) for every component of the amoeba complement, that is, for ev- 
ery monomial of the polynomial P(z,w) except for degenerate components, 
where the facet is reduced to a point. The slope of the facet corresponding 
to the term pijZ l w^ is always the same, namely The intercept, how- 

ever, varies. We will now show that these intercepts can be taken as local 
coordinates on the space of Harnack curves. 

Since the intercepts of the unbounded components are easily found from 
the coefficients on the boundary of the Newton polygon of P, and hence, 
from the boundary points of C, we can alternatively take the points on the 
boundary and the intercepts of the bounded facets as our coordinates. 

The space of Harnack curves with given boundary data and genus at 
most g is a real semialgebraic variety of dimension g, naturally stratified by 
the possible genus degenerations (where one or more ovals degenerates to a 
point). We have the following 

Proposition 6. The variety of Harnack curves with given boundary and 
genus is smooth with local coordinates given by the intercepts of the nontrivial 
compact ovals. 

Proof. The tangent space to a curve C with given boundary and nodes is 
formed by the space of polynomials R of degree d vanishing at all nodes and 
boundary points of C, modulo the polynomial P itself. Such polynomials 
have the form zw Q(z, w) where Q is a polynomial of degree < d — 3 vanishing 
at the nodes of C. It follows (see Section 12.2.4)1 that the space of such 
polynomials is ^-dimensional. 

Let Q be such polynomial and let (x, y) be a point inside the kth hole of 
the amoeba. We compute the variation of the R(x, y) in the direction Q as 
follows 



d_ 

dt 



Rp+tzwQ(x,y) 



\_ ff Q(z,w) 



dz dw 



(27rz) 2 JJ\z\=e x P(Z,W) 



w\=e : 



1 



/, - dZ E 

J \Z \=e X J-yj \ 




Q(z,w r ) 1 



Q(z,w) 



dz, (4.9) 



2ni 




P{z,w) 
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where Rp+tzwQ denotes the Ronkin function of the polynomial P+tzw Q, the 
integral over w is computed by residues, and fik is the 6-cycle corresponding 
to the /cth hole, see Section 12.2.41 

We see that the Jacobian of the transformation from the coefficients to 
the intercepts is the period matrix of the curve C, and, in particular, it is 
nondegenerate. □ 

Proposition |H1 can be strengthened as follows: 

Proposition 7. The neighborhood of genus g degree d Harnack curve inside 
the space of all degree d Harnack curves with the same boundary is diffeo- 
morphic to 



where the genus g stratum is embedded as (R 9 ,0). 

Proof. The R>o coordinates are given by the values of the polynomial Q 
from the proof of Proposition |U] at the nodes of C. These values have to be 
of definite sign in order for an oval to be present but otherwise arbitrary by 



4.3 Variational principle 

The holes of the amoeba satisfy the following variational principle which is 
parallel to the variational principle for conformal maps. 

Proposition 8. Under a boundary preserving variation that lowers one in- 
tercept while keeping all others fixed, the corresponding hole in the amoeba 
shrinks while all other components of the amoeba complement, including the 
unbounded ones, expand. 

The proof will use the Legendre transform R v of the Ronkin function. In 
the dimer problem, it has the meaning of the surface tension function. It is 
defined on the Newton polygon of P, satisfies the Monge- Ampere equation 



and has conical singularities (which correspond to facets of R) at the lattice 
points inside the Newton polygon corresponding to complementary compo- 
nents . 




Section l2~2~4l 



□ 




(4.10) 
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Proof. Let 5R V denote the variation of R v . It is positive at one singularity 
of R v and zero at all other singularities. It also vanishes on the boundary of 
the Newton polygon. 

The function 5R V satisfies the equation 

Rxx&Ryy ~t~ Ryy&Rxx Rxy^Ryx Ryx^Rxy 0; 

which is the linearization of (j4.10J) . and can be written 

ti(d 2 R y )-\d 2 5R w ) = U. (4.11) 

Here d 2 denotes the Hessian matrix of a function. Away from singularities, 
the function R v is strictly convex, so d 2 R v is positive definite, and hence 
()4.11|) satisfies a maximum principle: local extrema of SR V are not allowed. 
From the boundary conditions, we conclude that SR^ > everywhere. It 
follows that the cones at all vertices except one are becoming more acute, 
hence all but one of the facets of R are expanding. 

Adding a constant to i? v , which doesn't affect ()4.10|) . we can achieve that 
<5i? v is zero at one singularity and negative at all others. It follows that 
the cone at this one singularity is becoming more obtuse and, hence, the 
corresponding facet of R is shrinking. □ 

The following Figure El (best seen in color) illustrates the variational prin- 
ciple: 

4.4 Volume under the Ronkin function 

The maximum principle for the linearized Monge- Ampere equation used in 
the proof of Proposition |H] has another application for the volume under the 
Ronkin function. 

Let C\ and C% be two Harnack curves with the same boundary data and 
assume that their equations Pi and P2 are normalized in the same way (such 
as e.g. both constant terms are equal to 1). Then the corresponding Ronkin 
functions R\ and R2 agree asymptotically at infinity and the integral 

Voli?i - Vol i? 2 = // {Ri ~ R2) dxdy 
J Jr 2 

converges. Choosing any R 2 as our reference point, this defines a functional 
Vol-R of a Harnack curve. This function has the following monotonicity in 
the intercept coordinates. 
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Figure 5: One oval is shrinking, the others expanding. 



Proposition 9. A variation of the curve C which lowers the intercepts also 
lowers VolC. The genus zero Harnack curve is the unique volume minimizer 
with given boundary conditions. 

Proof. This follows from the same maximum principle used in the proof of 
Proposition |HJ a variation which lowers the intercepts necessarily lowers the 
entire amoeba. The intercepts can be lowered to the unique (by Theorem 
@]) point where the amoeba holes shrink to points, that is, to the genus-zero 
curve. □ 

Theorem 5. The variety of spectral data in (13.7)1 is diffeomorphic to M. d2+1 . 

Proof. The gradient flow of the volume functional with respect to any metric 
on the projective space of degree d curves contracts the variety of Harnack 
curve with given boundary to a small neighborhood of the genus zero curve. 
It follows that this space, together with its stratification by the genus is 
diffeomorphic to the product of ( n 1 ) copies of M> . Adding a standard 
divisor makes it diffeomorphic to ]R( rf-1 )( d ~ 2 ) (each hole, possibly degenerate, 
with a point on its boundary is diffeomorphic to an M. 2 ). Adding in the 
boundary points, which we showed in section 14.11 was a space diffeomorphic 
to R 3d ~\ makes R d2+1 . □ 
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4.5 Areas of holes 



Under the variation considered in Proposition |H] the area of one of the holes 
in the amoeba is decreasing, while the area of all other ones is increasing. 
Since the outer components of the amoeba complement are also expanding 
and the total area of the amoeba is preserved, we conclude that the decrease 
in the area of one hole dominates the total increase in the areas of all other 
holes. In other words, the Jacobian of the transformation from intercepts to 
the areas of holes in the amoeba is strictly diagonally dominant and, hence, 
nonsingular. We obtain the following 

Proposition 10. The areas on the amoeba holes can be takes as local coor- 
dinates on the manifold of Harnack curves with given boundary and genus. 

This can be sharpened as follows: 

Theorem 6. The areas of the amoeba holes map the set of all Harnack 
curves with given boundary diffeomorphically onto R> ^ 2 ^ 2 , mapping the 
stratification by the genus to the stratification by the number of nonzero co- 
ordinates. 

Proof. It remains to show that the map from the curve C with fixed boundary 
to the areas of its amoeba holes is proper. Suppose that some coefficients of 
the equation P(z, w) = J2ijPij z%w ^ °f C grow to infinity. Recall that fixing 
boundary means fixing the coefficients on the boundary <9A of the Newton 
polygon A, so only the interior coefficients Pij can grow. 

For each point (i, j) in the interior of A define the convex set 

Cij = {(a, b) eM 2 \ ai + bj + log Pij > max (ak + bl + \ogpki) }■ 

(k,l)^(i,j) 

The union over (i, j) of the C^j is the set 
{(a, b) G M 2 | max (ai + bj + logp.y) > max (ai + bj + logp^)}. (4.12) 

By our assumption, this set grows to infinity, in the sense that it eventually 
contains a ball of any arbitrarily large radius. Therefore one of the sets C^j 
eventually contains a ball of radius r, where r can be arbitrarily large. 

Since the Euclidean distance between (i,j) and any lattice point (k, I) ^ 
(i,j) is at least 1, for any (k, I) ^ (i,j) there exists (a , bo) satisfying a^ + bl < 
(r/2) 2 such that 

a {k-i)+b {l-j)>r/2. 
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It follows that if Cij contains a ball of radius r, the set of (a, b) satisfying 
ai + bj + hgpij > r/2 + max (ak + bl + \ogp k i) , (4-13) 

contains a ball of radius r/2. The inequality (j4.1Hj) implies that 

\pij z % w^\ > e r / 2 max \p H z k w l \ , \z\ = e a , \w\ = e b , 

therefore, all points (a, 6) satisfying (j4.13|) lie in the (i, j) component of the 
amoeba complement provided r is large enough. It follows that the area of 
that component is unbounded, which concludes the proof. □ 

Corollary 11. The areas of the amoeba holes and the distances between the 
amoeba tentacles provide global coordinates on the moduli space of Harnack 
curves. 

Note that the distances between the amoeba tentacles can be viewed as 
"renormalized" areas of the semibounded components of the amoeba com- 
plement. 

5 Isoradial dimers and genus zero curves 

Since Harnack curves of genus zero played a special role in our analysis, 
it is natural to ask for a characterization of those dimer weights that lead 
to spectral curves of genus zero. Genus zero weights are distinguished, for 
example, by maximizing the partition function for given boundary conditions. 

In this section, we show that, up to gauge equivalence, genus zero weights 
are the same as isoradial weights studied in [Hj. The results of [H] yield 
an explicit rational parametrization ()5.2|) of the spectral curve for isoradial 
dimers. In genus zero, the spectral curve determines the weights uniquely, 
up to gauge transformation. It remains, therefore, to characterize the curves 
(15. 2j) inside all genus zero Harnack curves, which is the content of Proposition 

ED 

5.1 Isoradial embeddings 

If T is embedded so that every face is inscribed in a circle of radius 1, we say 
that it is isoradial on condition that the weight of an edge is \/4 — £ 2 , where 
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Figure 6: Isoradial embedding of the honeycomb graph. 



£ is its length. Thus the weight is the distance between the circumcenters of 
the two faces adjacent to the edge. See Figure El 

For the hexagonal lattice, an isoradial embedding with d X d fundamental 
domain is determined by three <i-tuples of unit modulus complex numbers, 
{ax, . . . , ad}, {b\, ... , bd} and {ci, . . . , q}, which are edges of the rhombi con- 
necting vertices to the face centers as indicated in Figure El They obey the 
condition that they lie in three disjoint subintervals of the circle, that is, go- 
ing counterclockwise around the circle one encounters first the a's then the 6's 
then the c's. It is natural to consider these parameters up to a simultaneous 
rotation. 

5.2 Bloch-Floquet eigenfunctions 

In [TJ] it was shown that Bloch-Floquet eigenfunctions on isoradial graphs 
satisfy the "local multiplier condition", which says that if b and w are (re- 
spectively black and white) adjacent vertices on an edge bounding faces with 
circumcenters C\ and Ca, then 

[u - ri)(u - r 2 ) 
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where u is a complex parameter and = Ci — w are unit complex numbers. 
The constant y/r~ir~2 does not appear in ^1] due to the fact that we are using 
a different gauge for the Kasteleyn matrix. 

In particular the eigenfunction has Floquet multipliers 

/(w+ (1,0)) jju-Oi fb~ 

Z{U) = 7-- =11 T-\l ( 5 - 1 ) 

, v /( w + (0,1)) T-rW-Q /6~ 

/ ( w ) f = 7 u-bi\ ^ 

The parameter uGF 1 parametrizes the spectral curve P(z,w) = which is, 
therefore, rational. Setting 

u = e lt , Oi = e JQ!l , fei = e lft , q = e n " , 

the parametrization (j5.1|) becomes the following 

*)=n^. «(«)= ( 5 - 2 ) 

sm il ir 1 sm 

i 2 i 2 

in particular, t e R parametrizes the unique nontrivial oval. 



5.3 Characterization of isoradial curves 

Compare (|5.2|) to (|4.1j) and note the absence of arbitrary constant factors in 
front. We will call curves of the form (|5.2|) the isoradial curves. They have 
the following simple characterization 

Proposition 12. A genus zero Harnack curve C is isoradial if and only if 
the origin is in the amoeba of C . 

Proof. Putting u = in (|5.1|) we obtain a point on the unit torus, so the 
origin is in the amoeba of any such curve. 

Conversely, making the standard change of variable from the upper half 
plane to unit disk in Proposition we see that any Harnack curve of genus 
zero will have a parametrization of almost the same form as (j5.1|) . except 
we have to allow two arbitrary real factors in front. Such parametrization 
is unique up to the SX 2 (lR)-action by automorphism of the unit disk, so we 
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need to check whether the SL 2 (M.) -action can be used to set these arbitrary 
factors to 1. Conversely, we can act by SL 2 (R) on (|5.1|) and see what kind 
of constant factors we can get. 

Since rotating the unit disk around the origin clearly does not change 
anything, we have to look at the transformation 

u — C 

T(u) = -^, |C|<1. 
1 — Qu 

Applying it to (j5.1|) . we obtain the parametrization: 

z{u) = [I 



W(U) 



where a« = T 1 (aj) and similarly for 6, and 5j. Comparing this with (j5.1j) . 
we observe that 

and hence, by its very definition, the amoeba of the curve C describes the 
possible values of these factors. □ 

Note from the proof that the 2-to-l property implies that point (, when 
it exists, is uniquely defined. 

By a simple rescaling of z and w we can shift the amoeba of any curve 
C so that it contains the zero. In particular, for any genus-0 Harnack curve 
this gives a canonical family of isoradial parameterizations indexed by the 
points in the amoeba. 
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